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EVIDENCE FOB STRONGLY DAMPED GRAVITY 


WAVES IN THE EARTH'S ATMOSPHERE 

Igor J. Eberstein* 

Aeronomy Branch 

ABSTRACT 

Data from a series of pitot tube rocket soundings of the atmosphere are 
examined. The experiments consist of four soundings carried out over a two- 
day period at Fort Churchill, Canada in January and February 1967 by Smith, 
Theon and Horvath. 

Temperature data from the four soundings were averaged, and a smooth 
curve was drr.wn through the points. A hydrostatically determined atmosphere 
based on the above lapse rate v/as calculated. Deviations from the above mean 
atmosphere were calculated, and a wavelike structure was observed. The waves 
grow very slowly below approximately 80 km, and acquire a more rapid growth 
rate at higher altitudes. The density and temperature variations have a wave- 
length of 10-20 km and appear to consistently be about 180 degrees out of phase. 
The pressure variation is somewhat irregular. However, in regions where 
pressure variation is w^ell behaved it is out of phase with temperature and den- 
sity by about 90 degrees. 

Volland’s theory of thermally damped gravity waves was used to compute 
wave patterns like those observed experimentally. A highly damped, upward 
propagating wave gives good results. Damping by conduction, turbulence and 
radiative transfer is examined, and turbulence is found to be the most important. 


*NRC-NASA Research Associate 
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EVIDENCE FOR STRONGLY DAMPED GRAVITY 


WAVES IN THE EARTH’S ATMOSPHERE 
I. EXPERIMENTAL EVIDENCE 

Data from a series of pitot tube rocket experiments were studied. The ex- 
periments consisted of four soundings carried out over a two-day period at Fort 
Churchill, Canada in January and February 1967. The data are published by 
Smith et. al. (1967). The pitot tube technique is described by Horvath et. al. 
(1962). 

The above experiments measure stagnation pressure with a pitot tube. The 
velocity and position of the vehicle are also determined. Density and tempera- 
ture may then be determined by using the hydrostatic equation for the atmosphere 
and the ideal gas law. The underlying assumption of such a procedure is that 
vertical accelerations of air are small compared to the acceleration due to 
gravity. Present knowledge suggests that such an assumption is eminently 
reasonable. 

The temperature data from the four soundings were averaged and a smooth 
curve drawn through the points, as shown in Figure 1. Admittedly, some arbi- 
trariness is involved in drawing the mean temperature curve. However, such 
arbitrariness can only be satisfactorily resolved by reducing the scatter in the 
data. Such a reduction may be obtained by haring more rocket shots in a series. 

A hydrostatically determined atmosphere based on the above lapse rate was 
calculated. Experimental data were then plotted as p/p, P/P, T -T versus 
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altitude where the barred quantities refer to the mean atmosphere calculated 
above. The four plots, figures 2-5, exhibit consistent behavior. Temperature 
and density appear to consistently be about 180 degrees out of phase. The den- 
sity and temperature variations have a wavelength of 10-20 km. The pressure 
variation seems somewhat irregular. However, in the regions where pressure 
variation is well behaved it is out of phase with temperature and density varia- 
tion by about 90 degrees. The waves grow very slowly below approximately 
80 km, and acquire a more rapid growth rate at higher altitudes. 

A correlation analysis was performed on the data below 80 km. Figure 6 
shows auto- correlation plots for each of the four data sets individually, and for 
the four sets combined into one string of data. If the data is assumed to have a 
sinusoidal variation then the curve has the value 0.7 at V8, suggesting a wave- 
length of 12 km to 22 km with the mean at 16 km. If the variations are assumed 
to be random, then an integral scale of approximately 3.5 km is obtained. A 
power spectrum analysis of the combined density data was carried out, and is 
shown in Figure 7. The horizontal scale is wavelength and the vertical scale 
represents the spectral power density in arbitrary units. It is seen that wave- 
lengths between 2 km and 16 km make a significant contribution, with a peak at 
14 km. The analysis is not carried beyond 20 km wavelength since there is not 
enough data to justify such extension. For a discussion of correlation analysis 
and power spectra see Blackman and Tukey (1958), and Bendat and Piersol (1966). 
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II. GRAVITY WAVE THEORY 


Waves in a stratified fluid under the influence of gravity appear to have 
been initially discussed by Burnside (1889) and Love (1891). Both authors treated 
an incompressible fluid. Gortler (1943) used schiierin photography to show ex- 
perimentally that disturbances in an incompressible stratified medium under the 
influence of gravity propagate along characteristic rays. Lamb (1908) treated a 
compressible, adiabatic, ideal gas whose density is stratified by gravity. Since 
the original contributions by Burnside, Love and Lamb there have been many 
theoretical papers on various types of gravity waves. Reasonably up-to-date 
treatments of the subject may be found in Eckart (1960) and Yih (1965). How- 
ever, gravity wave theory is currently undergoing an active phase of develop- 
ment and contributions are being added to the literature at a rapid rate. 

The computer experiments tc be discussed in this paper are based on a 
gravity wave theory developed by Volland (1969). Volland assumes an ideal gas 
which is conductive but inviscid. Volland had previously concluded (Volland, 

1969, II) that, at least in the thermosphere, conductivity is more important as a 
damping mechanism than viscosity. In speaking of a conductive gas we shall 
mean a gas having an ’’effective conductivity,” i.e.,one in which a combination 
of thermal transfer processes behave in a manner which can be approximated 
by the heat conduction equation. Processes contributing to the effective con- 
ductivity include conduction, turbulent mixing, radiative transfer, and possibly 
chemical effects. The importance of these individual processes will be discussed 
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below. 


Since ‘tmospheric gradients are usually gentle* the ray approximation 
for vertical wave propagation was used. Volland (1969, 1) has shown that the 
ray treatment is a sufficient approximation for gravity waves propagating 
obliquely upwards, provided the horizontal wave number is in the order of, or 
greater than o/a where a is the sound speed, and cois the circular frequency of 
the wave. In the present study, the above condition was generally satisfied. 
The frequency of a moving fluid particle, or the ’intrinsic frequency’ may be 
defined as 


n - u> - k • u 


where 

k is wave number 
U is wind velocity 

is wave frequency in a quiescent atmosphere 

The above value of intrinsic frequency may be substituted into the equa- 
tions describing gravity wave propagation in a quiescent atmosphere, and a 
solution obtained. Jones (1969) discusses the conditions under which the above 
procedure is valid. Basically, the following assumptions must be satisfied: 

1. Vertical and horizontal derivatives of vertical winds must be less than 
(N 2 - fi 2 ) l/2 , where N is the Brunt- Vais ala frequency. 

2. Horizontal derivatives of horizontal winis must be much smaller than 0. 
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3. Time derivatives of horizontal winds must be much smaller than 
* (N 2 - w 2 ) l 2 , where ; : 1 N 2 g. 

4. Time derivatives of vertical winds are much smaller than g, the accelera- 
tion due to gravity. 

These conditions are generally well met, except when ft approaches either 
N or zero. 

III. ESTIMATES OF EFFECTIVE CONDUCTIVITY 

Processes contributing to "effective conductivity" include conduction, con- 
vection and turbulent mixing, radiative transfer and, possibly, chemical effects. 
The importance of these processes will now be considered. Conduction will be 
considered first. 

Define: 

?C p P 
" 2k 


* 


where 

y - ratio of specific heats 
C p = specific heat capacity of air 
k = conductivity of air 
Further, define 


G 
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where < is the circular frequency of the gravity wave whose propagation is Ix?- 
ing considered. The above parameters are due to Volland (1969). 

Now, the thermal diffusivity, * , is defined as: 



Substituting into the equation for one obtains: 



The term in the second bracket is the formula for ihe Pquare of the sound 
speed in an ideal gas. Thus one may write: 

a 2 

“h = 2 * ’ 

Figure 8 shows a plot of co h versus altitude for the standard atmosphere. 

It is seen that o> h , and therefore G, depend on density and temperature. Density 
variations of ±50% are readily observed, and similar variations in the damping 
parameters may be expected. 

Though G decreases exponentially with increasing altitude, its value is 
everywhere quite large, compared with values of about 40,000 necessary to ex- 
plain the slow upward wave growth which is observed. It is thus quite obvious 
that pure conduction does not provide an adequate damping mechanism, at least 
not below 120 km altitude. 
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The effects of conduction are generally enhanced by a combination of con- 
vection and turbulence. 

Assume that we may define an effective turbulent thermal diffusivity K , 
such that 

/Re - Re cfl( \ 

V Re o ) 

where * is the molecular thermal diffusivity, and Re is Reynold's number. The 
Reynold’s number is generally defined as follows: 

LV 

Re - — 

where L is some suitable length dimension, V is a suitable velocity and v is the 
kinematic viscosity. Following general astronomical practice, further justified 
by Blamont and DeJager (1961), L will be taken as the atmospheric scale 
height. Blamont and DeJager used a turbulent velocity, estimated from their 
measurements, for V. In the present study it will be assumed that 

~ „ — 

V “ flV 



where V the turbulent velocity, V is the mean wind velocity, and /3 is the ratio 
of r.m.s. turbulent fluctuation velocity to mean velocity of the flow, and is gen- 
erally in the order of 10%, All that is necessary for the following discussion is 
that (1 remain approximately constant. 
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Bl&mont and DeJager (1961) found that the atmosphere undergoes a turbulent 
to laminar transition at about 100 km altitude. At this point, the Reynold's 
number based on fluctuation velocity is 2,000. 

In the present study a mean wind structure having the following character- 
istics Is assumed: 


V = Yz 0 < Z < 5 

\ Z Z > 5 

v/here V is mean wind speed in meters per second, and Z is altitude in kilo- 
meters. The mean velocity is then used to calculate the Reynold’s number, 
resulting in a transition. Re of 20,000. This result agrees with Blamont, and 

'V — 

DeJager' s value if one substitutes V ~ 0. 1 V. 

Experimental data presented in the Meteorological Rocket Network Data 
Reports (1967), and the upper atmosphere wind measurements of Bedinger (1968) 
show that the assumed wind profile is reasonable. However, it must be realized 
that winds have great temporal variability, and therefore, any stationary model 
must be accepted with great caution. The wind variability of course suggests 
great variability in the turbulence level, and therefore, variations in the damping 
factor. 

The reference value K 0 was taken from early Eiffel tower measurements 
reported by Brunt (1952) for an altitude of some 150 meters. There, K 0 = 10 5 
cm 2 /sec. Since * 0 is 0.2 cm 2 /sec, the ratio, (K A) 0 = 5 x 10 s . 
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Figure 9 shows K * versus altitude. The two triangles iepresent K l /D 
where K l is turbulent mass diffusivity, and D is molecular diffusivity. The 
mass transfer data are from the work of Booker (1956) and Zimmerman and 
Champion (1963). The agreement is seen to be reasonable. Zimmerman and 
Champion's results also show approach to molecular behavior around 100 km 
altitude, in agreement with the results of Blamont and DeJager. 

Figure 8 shows *^ h based on turbulent heat transport versus altitude. Tur- 
bulent transport of chemically active species may further enhance the effective 
conductivity. 

Another contributor to effective conductivity is radiative transfer. Brunt 
(1952) showed than when photon mean path is short compared to a representa- 
tive length, the radiative energy transport behaves in a manner similar to con- 
duction, and may be treated by means of a conductivity equation. Brunt's con- 
clusion is what one would expect from the physical nature at the two processes, 
i.e., energy transfer by a molecular random walk is quite similar to energy 
transfer by a photon random walk. In both cases the mean free path must be 
small compared to a relevant physical dimension. For the region where water 
vapor absorption is important Brunt estimated a radiative diffusivity coef- 
ficient K p - 1.3 x 10 3 cm 2 /sec. The radiative diffusivity based on carbon di- 
oxide absorption will, of course, be smaller. 

Let us examine some of the properties of radiative diffusivity based on 
carbon dioxide. Elsasser and Culbertson (1960) review experimental data on 
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C0 2 absorption and conclude that below one atmosphere total pressure the low 
resolution absorption is relatively independent of total gas pressure, depending 
only on the amount of absorber gas. The mixing ratio of C0 2 in the atmosphere is 
approximately constant. The radiative conductivity is 


X Ih 

r 


where ■*. is mean absorption coefficient, and 


f a o * 


The radiative t' mal diffusivity is defined as 

X r 

K * ■ ^ C P ' 


It follows that 

P% _ *o 

K * = PC p = C p 

and is approximately constant throughout the atmosphere, provided the radia- 
tive properties themselves do not significantly change. 

Radiative transfer due to water and ozone are, of course, altitude dependent. 
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Golitsyn (1903) studied radiative damping of waves for long photon mean 
path. Golitsyn developed a coefficient of radiative conductivity given by: 



where ; is wave number of the radiation, and k is wave number of the gasdy- 
namic disturbance. % is the frequency dependent absorption coefficient, and B 
is the Planck function. 

Damping by carbon dioxide and ozone were studied, using low resolution 
absorption spectra given by Elsasser and Culbertson (1960). The results are 
shown in f igure 8. It is seen that for the longer vertical wavelengths radiation 
damping does become significant. However, the damping estimate in the band 
center is likely to be somewhat high for the longer wavelengths since there the 
mean photon path is not really long compared to the wavelength of the gravity 
wave. When the mean photon path is not long, then Golitsyn's formula tends to 
overestimate the attenuation. However, the mean photon path for shorter gravity 
waves as well as in the wings of the band for all wavelengths is long and a cor- 
rect damping estimate may be expected. Atmospheric ozone distributions were 
taken from Mitra (1952) and Krueger (1969). 

In conclusion it may be said that turbulent heat transport is the most im- 
portant damping mechanism with radiation making a significant contribution 
below 75 km. Above 110 km, conduction provides most of the damping. 
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IV. NUMERICA L CALCULATIONS 

The numerical calculations were based on the gravity wave treatment de- 
veloped by Volland (1909). Essentially, Volland uses a matrix method to obtain 
periodic solutions of a linearized version of the following set of equations: 

<> ■ j . , 

— i V • ( , V ) * 0 

■ 15T + VP * .* g " o 

DT 

/ C v op ^ P v * V - V(k VT ) 0 

P - mpT = 0 

Inclusion of the heat transfer term in the energy equation provides a 
mechanism of dissipation since destruction of thermal gradients leads to the 
creation of entropy. 

Solution of the linearized equations yields a dispersion relation which may 
be solved for the vertical wave number giving: 



where 

D = 


k 2 - io) 
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u ~ circular frequency of the wave 
a local sound speed 


k 


u 


) 


J B 

°A 

g 

H 

A 

B 

S 

G 


a 

ratio of specHic heats 
g a i , ~ \ ~ Brunt - Vaisala frequency 
/g 2a a 2H 
acceleration due to gravity 
Pressure scale height 


cc.- h - a* /2 k 

k - thermal diffusivity 

The adiabatic case is obtained in the limit as G becomes infinite. In that 
situation the dispersion relation becomes the same as obtained by Hines (1960), 


i. e. 



where 


n 


k a 

X 

co 


n 


, a 

t 

co 
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Georges (19(37) lias shown that the dispersion equation may be rearranged into 


the form 


k 


2 

x 


, 2 


t 


k 2 

l 



1 


where 


t 2 


2 

(J 




and 


m 


2 



As Co- -* GO 


-£ 2 


m 


2 


cu 


2 


a 


2 


giving: 



1 
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Ui 

a 



+ k f a 



k 


i.e. t the two-dimensional equation for a pure sound wave. 
Conversely, as w - o 



giving 




1 

4H 2 


* 


where use was made of the equality 


CO 

a 



_ 1 _ 

2H 


* 

The low frequency limit dispersion equation is very close to what would be 
obtained by using an incompressible, density stratified atmosphere with gravity 
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Tolstoy (1963) cited by Georges (1967), gives; 



For small ^ the Tolstoy formula becomes identical with the low frequency 
dispersion relation derived above. 

For 00 > a; > u> A the regime may be viewed as one of acoustic waves modi- 
fied by the influence of gravity. Clearly the gravitational influence will be most 
strongly felt for low frequency "infrasonic" waves. 

For a > a) > 0 the regime may be viewed as one of gravity waves modified 
by the effects of compressibility. 

Continuing the above form of resoning, the Volland dispersion relation rep- 
resents a compressible gravity wave modified by heat transfer. 

Originally a plane wave solution was assumed, and shown to be compatible 
with the linearized equations. One may solve for the physical parameters to ob- 
tain what Hines (1960) has called polarization relations. The physical parameters 
may then be written in matrix form as 

(e) = (K) (c) 


where 


is the matrix of physical parameters with 


where 




u_ S b 

a ~ e a y 


AT 

T = e, i 


b 


i 



is an altitude normalization factor. 

w is vertical perturbation velocity 
u is horizontal pertubation velocity 
c is the "basics wave" 

c = c 0 exp j^i (cot - k f * - k x x 
Taking the wave at t * 0, x = 0, one obtains 

k, s 

c = c 0 e * 

for an isothermal atmosphere. 
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The ray solution for a series of isothermal slabs then gives 


where 


c 


c 0 exp 



Q(C) d£ . 


The transformation, or polarization, matrix K is given byi 


where 



Lj = q - iA«< 


1 - 


2S 2 


L, = y - 2A 2 - 2i Aq 


L = 7 - S 2 - A 2 - q 2 


The above matrix of Volland (1968) reduces to Hines* (1960) polarization 
relations for ' 3 adiabatic case. 

A mean temperature profile was fed into the computer program which then 
calculated a mean atmosphere and various characteristic quantities based on the 
mean atmosphere at 1 km vertical intervals. The atmosphere was assumed to 
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be homogeneous horizontally. The mean atmosphere was based on the hydro- 
static equation and the ideal gas law. 

Values of the circular frequency, < , and the horizontal wavelength were 
selected. The range of frequencies considered was 10" 3 sec to 10" 2 sec. The 
horizontal wavelength, K n was varied between 25 km and 500 km. 

The initial amplitude of the elementary wave was adjusted to match the 
density deviation in the 20 km to 40 km altitude region. 

A series of values of the damping parameter G were selected. A fairly 
high value of damping was found to be necessary to give a wave growth compar- 
able to what was observed experimentally. 

Most calculations were for a quiescent atmosphere, but several wind dis- 
tributions were also studied. A wind which varied sinusoidally with altitude 
was found to give the most interesting results. 

V. DISCUSSION OF RESULTS 

Figure 10 shows density fluctuation for a series of gravity waves, the first 
two of which have damping characteristics shown in Figure 8. The third curve 
has the same variation of damping coefficient with altitude, only damping is 
twice as large as in the previous two cases. Figure 11 shows a situation where 
the value of G is kept constant at 40,000. For the value of omega chosen 
(9 x i , " 4 sec" 1 ) the above implies w H - 36 sec" 1 » constant, i. e. approximately 
twice as much damping as estimated in Figure 8. Figure 12 show? the situation 
where the damping coefficient increases linearly with altitude, and G - G 0 /Z 
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where Z is altitude in km, G 0 - 940,000. For u> = 9.0 x 10” 4 see** 1 , = 18 sec” 1 

at 50 km and 8.5 sec* 1 at 100 km. These values imply an order of magnitude 
more damping than estimated. Whereas the damping estimate could be too low, 
it i» more likely to be on the high side. Thus, different explanations are called 
for. It is also unlikely that contributions from viscosity effects would be an 
order of magnitude more significant than the thermal effects. Two alternatives 
remain. Either the damping is due to non-linear terms not considered in the 
theory, or additional turbulence is generated by the gravity waves, so that 
eventually the waves are propagating into a turbulence field which is much 
stronger than the mean field existing in the absence of gravity wave activity. 
Such a conclusion is eminently reasonable, considering that a gravity wave 
having a 30% density deviation generates a wind well over 100 m/sec. 

Gossard (1962) observed gravity waves in the troposphere. He found the 
perio ! to vary between 15 minutes and 2 hours. The horizontal wavelength fer 
waves of 15 minute period was 19 km and for waves of 2 hours period it was 
150 km. The periods investigated in the present study range from 19.5 minutes 
to 1.94 hours, thus being approximately the same as those observed by Gossard. 
However, the horizontal wavelength range was between 372.5 km and 48.5 km, 
i.e., about 2. 5 times the wavelengths reported by Gossard. The difference is 
due to the longer vertical wavelengths which are indicated by the present study. 
Quite approximately 
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Thus, for a givon atmosphere and period the vertical wavelength is approx- 
imately proportional to the horizontal wavelength, and a larger vertical wave- 
length calls for a larger horizontal wavelength. 

The sinusoidal background wind field shown in Figure 12 results in con- 
siderable modulation of the gravity wave. The modulation is of the same type 
as observed in the experiments, only more pronounced. 

The winds correspond to typical conditions both as regards amplitude and 
wavelength. Figure 13 shows an experimental wind pattern obtained at Fort 
Churchill by Smith et al. (1965). 

CONCLUSIONS 

Semi-quantitative agreement between an atmospheric wavelike structure 
and gravity wave theory has been demonstrated. Also, the study indicates what 
further experiments will be needed to close the loopholes in our understanding. 

The most important damping mechanism is turbulence. However, naturally 
occurring turbulence does not appear adequate to provide the experimentally re- 
quired damping It is suggested that gravity wave activity creates its own, much 
stronger turbulence field. Thus it would be important to measure atmospheric 
turbulence both at times of low and high gravity wave activity. 

Modulation of the wave by the wind Held is important. Thus, it would be 
important to simultaneously measure the thermodynamic parameters and the 
background windfield. 
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APPENDIX A 


The following contains some further thoughts on turbulent damping. In 
Figure 14, samping estimates by various authors are shown. 

Hosier (1969) determined vertical diffusivity from radon profiles. Average 
derived vertical diffusivities for a 90 meter height interval range from about 
3.1 x 10 5 cm 2 /sec during unstable convective periods to about 8.3 x 10 3 cm 2 /sec 
during stable inversion periods. The above measurements were carried out at a 
tower site near Washington, D. C. The Eiffel tower measurements reported by 
Brunt, and discussed in the main text, gave 10 5 cm 2 /see which is close to Hosier's 
value for unstable convective periods. 

The estimates of Johnson and Wilkins (1965, 1966) are based on estimated 
limits of steady downward heat transport in the atmosphere. Colegrove, Johnson 
and Hanson (1966) estimate an average eddy diffusion coefficient based on molecu- 
lar and atomic oxygen transport. 

The experimental damping curve shows the amount of damping necessary to 
maintain the observed growth rate of gravity waves with altitude. It is seen that 
the needed damping is considerably greater than any steady state estimates. 
However, the experimental damping estimates are below the diffusion coefficients 
determined by Hodges (1969) for a gravity wave where amplitude does not increase 
with altitude. 
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Figure 1. Mean Temperature Fro ,i Four Rocket Soundings 
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Figure 2. Deviation of Atmospheric Density, Pressure, and Temperature From The Mean 
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Figure 4. Deviation of Atmospheric Density, Pressure, and Temperature From The Mean 
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Figure 5. Deviation of Atmospheric Density, Pressure, nno Temperature From The Mean 
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Figure 9. Estimated Ratio Of Turbulent To Molecular Diffusivity As A Function Of Altitude 


34 



THERMAL 

DAMPING PARAMETER 



30 40 50 60 70 80 90 100 110 120 


ALTITUDE (km) 


Figure 10. Damping Parameter, G # As A Function Of Altitude 
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Figure 11. Calculated Gravity Wave Propagation In A Quiescent Atmosphere 
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Figure 14. Estimates of Eddy Diffusivity 




